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Abstract 

We reanalyze the dissociation process of the J/ip by tt and p mesons into 
D + D, D* + D, D + D*, and D* + D* within a meson exchange model. 
In addition to the dissociation mechanisms considered in the literature, we 
consider anomalous parity interactions, whose couplings are constrained by 
heavy quark spin symmetry and phenomenology. This opens new dissocia- 
tion channels and adds new diagrams in the previously considered processes. 
Compared to the previous results, we find that these new additions have only 
a minor effect on the p + J /ip total inelastic cross section, but reduce the one 
for 7r + J /ip by about 50 % near the threshold. 
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I. INTRODUCTION 



Suppression of J/ip production is considered to be one of the promising signals for de- 
tecting the formation of the quark-gluon plasma (QGP) in relativistic heavy ion collisions 
(RHIC) . Indeed, recent reports by the NA50 Collaboration |||| show an anomalous 
suppression of J/ip production in Pb+Pb collisions at CERN. Model calculations, which 
assume the formation of a deconfined state of quarks and gluons []5]|| at the initial stages 
of the collision, seem to be successful in explaining the anomalous suppression. However, 
before coming to a definite conclusion on the existence of the QGP state in such collisions, it 
is essential to investigate whether the observed J / ip suppression can be explained by a more 
conventional approach that does not assume the existence of QGP state. In such models, 
the suppression comes from the absorption of J /if) by the co- moving hadrons |7|-[nj. Since 
such co-mover models strongly depend on the cross sections of J / if) absorption by hadrons, 
it is essential to have a better knowledge of the magnitude and energy dependence of the 
cross sections. This is not an easy task, but is indispensable for true identification of QGP 
in RHIC. 

The absorption of J/if> by light hadrons has been investigated using various methods 



TT|— pT8| . The estimated values for the cross sections show a strong dependence on the 
hadronic models and assumptions on the absorption mechanisms employed in the calculation. 
The first set of models uses the quark degrees of freedom. Using perturbative QCD in the 



heavy charm quark mass limit |T§|] , Kharzeev, Satz, and their collaborators [0,0 found 
very small absorption cross sections at the order of ph. In Ref. ||13|| , Martins, Blaschke, and 
Quack investigated the charmonium dissociation processes, tt + J /if) — > D* +D, D +D*, and 
D* + D* using quark exchange model and found a strong enhancement near the threshold 
with an exponential falloff at higher energies. The peak was found at yfs ~ 4 GeV with a 
value of about 7 mb. This calculation was recently improved by Wong, Swanson, and Barnes 
14| by using a more successful quark-interchange model and hadron wave functions. The 



cross section for n + J/ip was found to be relatively small with a maximum of about 1 mb, 
while the pion induced if)' dissociation process had a much larger cross section with lower 
threshold for the initial kinetic energy. 

Another approach to J / if) absorption is to use effective meson Lagrangian [|15|-f[8|| . In this 
approach, the authors considered t-channel meson exchanges to estimate the J / if) absorption 
cross sections. In Ref. ||15| , Matinyan and Miiller calculated D meson exchange diagrams 
and found small cross sections around 0.3 mb for rr + J/ip and p + J/ip. The obtained 
cross sections were found to increase with energy. This calculation was improved by Haglin 
16| , who included D* meson exchanges and four-point couplings in the J /if) absorption 



mechanisms. Similar calculations were recently performed by Lin and Ko W%. The new 
results show a rather large cross section for 7r + J/ if) — > D* + D and D + D*, i.e., about 
20 mb at y/s = 4 GeV, which at higher energies saturates to about 30 mb.Q However, for 
a realistic structure of the mesons, form factors have to be taken into account. Inclusion 



1 The typo errors in Ref. and different conventions between Refs. [1G . 17 1 are cleared recently 
by Ref. p§. 
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of form factors generally reduces the cross sections. With suitable form factors and cutoff 
parameters, which however cannot be justified a priori, Lin and Ko [17| concluded that the 
saturated absorption cross sections for it + J/ip and p + J/ip are around 7 mb and 3 mb, 
respectively. 

All of the meson exchange model calculations discussed above essentially use the same 
effective meson Lagrangian; namely the minimal SU(4) Yang- Mills Lagrangian. Two imme- 
diate questions can arise. The first is the use of SU(4) symmetry in constructing the effective 
Lagrangian. Obviously, charm quark is much heavier than the other three light quarks and 
the SU(4) symmetry is nowhere near being true in QCD. Nevertheless, the starting point 
of using the SU(4) Lagrangian is to categorize the possible interaction vertices among the 
meson multiplets and estimate as many of their respective couplings by phenomenology as 
possible. Doing so, not all of the couplings can be fixed. However, it turned out that for 
couplings that can be checked, one finds that the SU(4) symmetry relations for coupling 
constants are not totally meaningless |I7,IB||. Such checks provide the order of uncertainty 
for the calculated absorption cross sections within the effective Lagrangian method, and 
makes the calculation meaningful. The second question is whether all possible interactions 
are considered. It is to this question that we want to address in this paper. The previous 
studies did not fully take into account all possible set of interaction Lagrangian involving D 
and D* . In fact, the anomalous parity interactions which are connected to the gauged Wess- 
Zumino action are all missing. Furthermore, the anomalous term such as D*D*n is essential 
and required by heavy quark spin symmetry and its coupling is strongly constrained to the 
D*Dtt coupling |20|pT| . As we shall see, inclusion of D*D*n coupling opens new absorption 
channels and mechanisms which were not considered in the literature. Therefore, to fully 
estimate the J/ip absorption cross section within the hadronic model, it is necessary and 
important to see the effect of the anomalous interactions to such processes, which we will 
investigate in this work. 

This paper is organized as follows. In Sec. II, we introduce our effective Lagrangian. 
The interaction Lagrangian needed for J /if) absorption processes by 7r and p mesons are 
obtained and the cross sections are evaluated. Then in Sec. Ill, we determine the coupling 
constants and give the numerical results for the cross sections. Section IV contains summary 
and discussions. Some details on deriving our effective Lagrangian are given in Appendix. 



II. EFFECTIVE LAGRANGIAN FOR J/ip ABSORPTION 

In the heavy quark mass limit (m,Q — ► oo), effective Lagrangian for light (qq) and heavy 
(q~Q or Qq) mesons can be constructed to preserve chiral and heavy quark symmetries 
p0^p4[. When the heavy quark flavor is the charm, finite mass corrections are expected 
to be important and should be introduced in the effective Lagrangian in a systematic way. 
The couplings of the finite mass corrections can be fixed by "velocity reparameterization 
invariance" |25| or by phenomenology. This "top-down" approach seems to be most plausible 



way to construct the effective Lagrangian. However, it is not yet clear how to implement 
heavy quark symmetry to the interactions of quarkonia with light hadrons as the quarkonium 
states (QQ) such as the J/ip contain two heavy constituents. 



Another approach which was used in previous investigations |15-18| is the "bottom-up" 
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approach. In this approach, one starts with the SU(4) symmetric chiral Lagrangian. The 
heavy charm quark mass is then assumed to be taken into account through appropriate 
symmetry breaking terms | 26fl . However, since the SU(4) symmetry is badly broken, finding 
all the appropriate symmetry breaking terms is very hard and beyond the scope of this work. 
In this paper, therefore, we use the "bottom-up" approach phenomenologically to categorize 
all the possible interaction vertices among the meson multiplets. The couplings are then 
determined by experimental data and by implementing heavy quark spin symmetry. Doing 
so, most of the three point couplings can be determined. For the four point couplings and 
some three point couplings, where such determination is not possible, we will partly rely on 
the SU(4) relations as well as other model predictions. 

As in Refs. |l6Hl8|, we start with the SU(4) Lagrangian for the pseudoscalar mesons and 
introduce the vector and axial- vector mesons by the massive Yang-Mills approach . The 
axial- vector fields can be gauged-away as in Ref. [p8fl . The procedure to obtain the effective 
Lagrangian is well explained, e.g., in Refs. [0,0] and is summarized in Appendix. The 
effective Lagrangian for tt, p, D, and D* then reads 



£>D*Dn 
£tpD*D* 

£-ipD*Dn 
£>DDp 

Co*D*p 



ig^DD^ (d^DD - Dd»D) , 

9i,n*D„r {dttd; + D;ttD) , 
Wddp (Dp»d,D - d,Dp»D) , 



1 9d*d* p 



£ipDDp 
£ipD*D*p 



where ir 



d u .D* v p a D* v - Dlp^D* 
+ D*» (p»d,Dl-d,p u D™) 

<?^*d* A (2D* V P"D* U - D*»p v D*» - D^ffD*) , 
r • 7r, p = t ■ p, and we define the charm meson iso-doublets as 



(la) 
(lb) 

(lc) 

(Id) 
(le) 

(If) 

(lg) 
(Ih) 



D T = 



(D°,D-), D=(D Q ,D+), 
_ (£)*o D *-j ; D * = ( D *a D * 



(2) 



In addition to the normal terms given above which were used by Refs. |16H18|, there 



are anomalous parity terms. The anomalous parity interactions with vector fields can be 
described in terms of the gauged Wess-Zumino action [p9] and are summarized in Appendix. 
By the same method used to obtain the normal terms, we get 



£d*d*tt 
£>4>d*d 

£-ipDDn 



-g^D^d^v {d a Dp + Dd a D*p) 
-ig fDD ^ v ^d v Dd^dpD : 



(3a) 
(3b) 
(3c) 
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C^ D ^ = -ig^^e^^Dtd^D* - ih^^e^d^D^D*, (3d) 
£d*d p = -9 D * Dp e^ al3 {Dd,p u d a D; + d,D* u d a p p D) , (3e) 

~ ih^^e^, (Dp„d a D; - d v D* aPf3 D) , (3f ) 

with £0123 = +1- In the heavy mass limit, due to the heavy quark spin symmetry, the 
pseudoscalars and vectors become degenerate and have to be treated on the same footing. 
As a consequence, the anomalous D*D*tt interaction term is required and can be related to 
the D*Dir interaction term by the heavy quark spin symmetry [p0| , plj| . 

Let us first evaluate the absorption processes of the J/ip by it and p mesons within the 
effective Lagrangian obtained above. The values of the coupling constants will be discussed 
later. In this study, we assume that the J ftp is a pure cc state and consider the OZI- 
preserving processes only at the tree level. The diagrams we are calculating are shown 
in Figs. [I] and |[ Compared to the previous studies Hl6lJl7|, we find that the anomalous 
parity interactions open new absorption channels, namely ir + J/ip — ► D + D [Diagram (1)], 
7r+ J/ip — > D* + D* [Diagram (3)], and p+J/ip — ► D* + D [Diagram (5)]. These processes are 
not allowed by the normal terms and were not considered before in the effective Lagrangian 



approach. In the quark exchange model 1131,141) however, the processes ir+J/ijj^D*-\-D 



and p + J/ip — > D* + D were calculated. There, ir + J/ip — ► D + D was not considered since 
the absorption amplitude vanishes unless one considers the spin-orbit force |T4] , |30| . A recent 
study using the Dyson- Schwinger formalism shows that the cross section for ir+J/ip — > D+D 
is small []30|. The anomalous terms also allow new absorption mechanisms to the previously 



considered absorption channels. These are Diagrams (2b), (4c), (4d), (6a), and (6b) of Figs. 
I] and [| As we shall see, the role of those diagrams is not negligible especially for the J/ip 
absorption by the 7r's. 

We first calculate the absorption amplitude of it + J/ip. We define the four-momenta 
of the J/ip and 7r by p\ and P2, respectively. We also denote the four-momenta of the final 
particles by ^3 and P4, which then defines s = (pi +P2) 2 and t = (p± —ps) 2 - The polarization 
vector of the vector meson with momentum pi is represented by e t -. Md and Md* represent 
the D and D* meson masses, respectively. Then the differential cross section for this process 
is 

lb 1 D^i 2 . w 



2 



dt 967rspf 



spm 



where pi is the three-momentum of p\ in the center of mass frame. Note that we have mul- 
tiplied factor 2 to get Eq. (^). This comes from summing over the possible isospin quantum 
numbers of the final state jTTj. In computing Diagram (1) in Fig. [I], the four-momenta of 



the D and D are represented by p% and p$, respectively. The absorption amplitude for this 
diagram can be written as 

M^MpM, (5) 

i 

where 
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M { } a) 



9D*Dn9^,D*D f aj3 _ (P2-P3) a (P2-P3)' 3 



(ft. - Pi) 2 - Mi: E ^ Xp \ 9 Ml* f V2 " P 1 Pi 

9d*Dtt9^d*d j ap _ (pi -Pz) a {pi -PzY 



M ( J C) = g^DDnS^XpP^PsPl. (6) 

In Diagram (2) of Fig. |l], the four-momenta of the D* and D are denoted by p^ and p±, 
respectively, and we have 

X=^E#?- (7) 

i 

where 

..(26) 9d*D*k9^D*D / afl (P2 ~ Ps) a (P2 ~ Ps)" \ 7 p A <5 



i>(2c) = 1 9d*d^d*d* \ af3 (pi-PsTiPi-Psf ] 
V " 2( Pl -p3) 2 -M|„ r Mf,* J 

X (p 2 a + P4a) [2^3^/3 ~ (Pi + Ps)/^ + 2 Plv^] , 

M™ = -g i!D * D ^v. (8) 

For Diagram (3) of Fig. [I], we write the four-momenta of the D* and D* by p% and £>4, 
respectively, and we have 

M = e?e?J2 M ^x £ i> ( 9 ) 



where 



i A (3a) 9d*D-k9x!)D*D 7 5 

M ^ = (P,-P3?-Ml £ ^ f ^ 



..(36) 9d*D-k9tI>D*D 7 5 

. .(3c) _ _ 9D*D*n9^D*D* j af3 _ (>2 ~ Ps)"^ ~ P?,Y 



X e^aPg^ ^4^/3 A - (Pi + P4)/9^A + 2piA^] , 
M (3d) 9p*D*^D*D* ( a p _ (gi -p 3 ) Q (Pl -P3)^ 

W (Pi-P 3 ) 2 -M|„ \ 5 M^, 

X £A 7 faPM [ 2 P3p.9f3v - (Pi + Pa)/^ + IPtogyfl] , 
•M-fjtuX = ~9iPD*D*tt £ ^pP2 ~ ^D*D*-K E Hv\pP\- (10) 

We now calculate the J j ip absorption by the p meson as shown in Fig. 0. We first define 
the four- momenta of the J/ip and p as p\ and p<i respectively. Then the differential cross 
section reads 
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d 4= 1 2 y\M\\ (ii) 

dt 2887rspf^' 1 v ; 

spin 

For Diagram (4) of Fig. |2|, we define p 3 and p 4 as the four-momenta of D and D respectively. 
Then, 

M = J2M^e^ 2 , (12) 



where 



JW ( 4fl ) _ A 9 D Dp9^DD 
W(46) A 9DDp9j,DD 

"* ~ (ft - Ps) 2 - JWJ. r J^WPiftftP* 

^ (4d ) _ 9D'Dp9^D'D ( af3 _ (Pl-P3) a (p2-P3) (3 \ ^.^.7 

-^t = 9^dd p 9^- (13) 

For Diagram (5) of Fig. 0, p 3 and p 4 are defined to be the four-momenta of the D* and 
D respectively. We have 

M=e?J2 M %i £ i^ ( 14 ) 



where 



/U (5a) - 2 9d*d p 9^dd s 7 

~ (P2 - Ps) 2 - M|, " A7i 

(5b) _ ^9DDp9^D*D 7 <j 

^A- (pi _ p3) 2_ M 2^PlP 3 P4„ 



Ayf (5c) _ 9 D *D*p9^D*D j a/3 _ (P2 ~ PgHgg ~ Pa)^ 



(ft-P3) 2 -M^r mi,. 

X Zp-fSpPjpl [2p3v9\a ~ (P2 + P^aQvX + ^2X9, 



va\ i 



u (5d) _ 9 D *Dp9^D*D* j a/3 _ (Pi -Ps^Pl ~ Ps)^ 



X i^foP^I [ 2 P3/^/3A ~ (Pi + Ps)/?^ + ^Pl\9^] , 

■Mpvx = ~9ipD*Dp £ vv\sp\ + h^ D * Dp ep U \sp s 3 . (15) 

Finally for Diagram (6) of Fig. |2|, we denote the four-momenta of the D* and D* by p% 
and p 4 , respectively, then we have 

M = eTef^2 M tX £ i^ ( 16 ) 
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and 



kA&o) 9d*d p 9^d*d a s 7 p 

JVi y,v\p — ^ —p 3 ) 2 — M 2 £yX ~< s£ w^P^^Pa- 



A/f (6b) _ 9 D *Dp9^D*D rflr^rfirfi 
JVL pu\p — ^ — p 3 ) 2 — Mp t ^<*(lP\P2P?,P^ 

M (fSc) 9 D *D*p9^D*D* f af3 (p2-P3) a (P2-P3Y 

M-Ap {p2 _ p3) 2_ M 2 f ]9 ^ 

X [2p 3 ug a \ - (P2 + P3)a9u\ + 2p2\9ua] \^P^9f3p ~ (Pi + Pi)(39pp + 2^1^^/?] 

w(6d) 9d*d* p 9^d*d* ( a/3 {pi-pzTipi-pzf 

W " (Pi - Ps) 2 - M 2 D « \ 9 M 2 D , 

X [2p 4l /#ap - (p 2 + P4)a0«/p + 2p 2 p5w] [2p 3 p5'/3A - (Pi + P3 )/30/*A + 2pVi9tf] , 

•^pJap = -9i,D*D* P {^g^gxp - g w g v x - g»\g vp ) ■ (17) 

In writing the above amplitudes, we have used pi ■ e$ = 0. Note that the cross section for 
7r(p) + J/^ ^ D + D* is the same as that of vr(p) + J/^ -> D* + £>. 



III. COUPLING CONSTANTS AND CROSS SECTIONS 

A. Coupling constants 

To estimate the cross sections, let us first determine the coupling constants of our effective 
Lagrangian ([[]) and (0). We follow the methods of Refs. fl6| , |r7f to determine the couplings 
for the normal interactions. Here we briefly explain the method referring to Refs. [|l6l , |T7i for 
details. The coupling constant of D*Dtt can be determined from the experimental data of 
D* — > Dn decay. Our effective Lagrangian gives 



D*^Di\ 



g 2 D*D-K I Ptt 1 3 
24vr M 2 d *' 



(18) 



However, at present, only its upper bound (~ 0.131 MeV) is known experimentally |3ll , |32| , 
which translates into g^*^ < 14.7. Therefore, we can only rely on model predictions such 
as those based on the QCD sum rule approach |33[] or relativistic quark potential model |34 



Here we will take the value obtained in Ref. |33[ and useQ 

9d*d« = 8-8- (19) 

For i[)DD, DDp and tpD*D*, D*D*p couplings, we will follow Refs. |l5| , p7|| and make 
use of the vector meson dominance model (VDM). By applying VDM to the couplings of 
r )D + D~ and / -fD°D°, one obtains fll~7 



2 Note that our convention is different from that of Ref. [33] by a factor of v2 
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2 1 

9tpDD = 7^9 J/4, = 7 - 71 ; 9dd p = -j9 P = 2 - 52 ' ( 20 ) 

where g v is determined from 

_ 47ra e 2 m M v 

1 y^ e + e - r , [ZL) 

° 9v 

with the vector meson mass My and a em = e 2 /4ir. The same method can be applied to the 
couplings of ^D* + D*~ and ■yD*°D*°. This gives 

9^D*D* = QijjDDi 9D*D*p = 9DDp- ( 22 ) 

It is shown in Ref. [17] that the SU(4) symmetry values are not far from the above estimates 
except g^ DD . 

There is no experimental or phenomenological informations on the 4-point vertices. 
Hence, we will rely on the SU(4) symmetry relations. Following Ref. |l7j, therefore, we 
assume that the appropriate symmetry breaking effects in the 4-point coupling constants 
are encoded via their relations to the 3-point vertices within SU(4) symmetry. Hence, using 
the phenomenological estimate of the 3-point vertices fixed above, we have 

1 

9i(,d*Dtt — ~^94,dd9d*d-k ~ 33.92, 
9ipDDp = ^94>dd9ddp ~ 38.86, 
94>d*d*p = 94,d*d*9d*d*p ~ 19.43. (23) 

We determine the couplings of the anomalous interactions in a similar way. As discussed 
before, the D*D*tt coupling is constrained by the heavy quark spin symmetry. Comparing 
with the effective Lagrangian of Ref. ||24j| , we get to leading order in 1/M(D) 

9d*d*. = ^p-9 D *D* « 9.08 GeV" 1 , (24) 

where M(D) is the average mass of D and D*. 

For ipD*D and D*Dp couplings, we can apply VDM to the radiative decays of D* into 
D, i.e., D* — > Z>y. Then by using the same technique explained before, we obtain 

2 ± _ 2 o 
9i,D*D — ~^9j/i,9v ~ ~^9j/i>9v> 

9d-dp = 19 p {9v + 2 9v)> ( 25 ) 



where g v '° are defined through 



(D+(k)\j; m \D* + (p,e)) = eg±e, uaP e v ky^ 

(D {k)\J° m \D*°{p,e)) = eg v e, ua pe»k a pP. (26) 

Then from the experimental values for the D* radiative decays, we can determine the above 
coupling constants. However, at present, only an experimental upper bound is known for 



9 



the D* decays. Therefore, we will use model prediction based on relativistic potential model 
given in Ref. [34], which gives 

* = £ + £ (27) 

where e Q is the heavy quark charge and e q the light quark charge with Aq = 1.57 and 
A g = 0.48 in GeV unit. This gives 

0^ « 1.12 GeV~\ $ « 0.96 GeV -1 , (28) 

which then leads to 

9ipD* d = 7-40 ~ 8.64 GeV" 1 , g D , Dp = 2.82 GeV- 1 . (29) 
For the 4-point couplings, we depend on the SU(4) relations again, which gives 

9tpr)* r)*n - h^ D * D * n = ^9^d*d9d*d-k ~ 38.19 GeV -1 , 

9ipD*D P = h"ipD*D P = 9^d*d9d*d* p ~ 21.77 GeV . (30) 

However, for the ipDDrr coupling, it is not easy to write it in terms of the other 3-point 
coupling constants as can be seen from the difference in dimensions. Thus we directly use 
the SU(4) relation and assume that the symmetry breaking effects change F n to Fjj |3l 



9tpDDiT = ^f^T « 16 - GeV_3 - ( 31 ) 



Here we have used Fd ~ 2.3F„- |[32|| . This completes the determination of our coupling 
constants. The values of the coupling constants are listed in Table | together with their 
fixing procedure. 



B. Cross sections without form factors 

We first give the cross sections for the J / ip absorption processes without any form factors. 
Since our meson exchange model is not expected to work at high energies, we focus on the 
energy region near threshold in this paper. In Fig. |3| we show the cross sections for the J / ip 
absorption processes obtained from Figs. [I] and 0. Although the ix + J/ip —>■ D + D process 
has the lowest threshold energy, its cross section (dotted line in the upper panel of Fig. ^) 
is suppressed compared to the other channels over the whole energy region. Nevertheless, 
this process was not considered previously in the meson exchange models. When the energy 
is less than the D* + D threshold, this process alone is allowed, but the magnitude is very 
small and less than 0.1 mb. This observation is consistent with Ref. [SC], although its energy 
dependence is different. 

The 7r + J /if) — ► D + D* has the same cross section as the tt + J/iI)^D* + D, and their 
sum is given by the dashed line in the upper panel of Fig. |^. These processes dominate 
for the energies from their threshold up to the D* + D* threshold. At higher energies, the 
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n + J/tp — > D* + D* process dominates the absorption. Since the D* + D* threshold is larger 
than 4 GeV, the n + J/ip^D* + D,D + D* is the dominant process for yjs < 4 GeV. 

The cross section for 7r + J /if) — ► £>* + D* process (the dot-dashed line in the upper panel 
of Fig. is the dominant one at large energies, which is in contrast to the quark exchange 



model calculations ||13|J14}| . In addition to the effects from the form factors, this may indicate 
that the meson exchange approach for this reaction has non-trivial contributions from the 
exchanges of meson resonance with heavier mass, such as the Di exchange, because of its 
higher threshold. This process is allowed by the anomalous terms and was not considered 
in the previous meson exchange model studies. 

In the lower panel of Fig. ^ we show the cross sections for the p + J /if) processes. 
It is shown that at very low energy, the exothermic p + J/ if) — > D + D is the dominant 
process (the dotted line in the lower panel of Fig. ||). But it is rapidly taken over by 
p + J /if) — > D* + D,D + D* process (the dashed line in the lower panel of Fig. |||), which 
also dominates at higher energies. This reaction is an anomalous parity process and was not 
studied before in the meson exchange model. The p + J /if) — ► D* + D* cross section (the 
dot-dashed line in the lower panel of Fig. has a peak near threshold but has comparable 
size with that of p + J/ip -> D* + D,D + D*. 

The dashed lines in Fig. f| show the results obtained with the normal terms only and 



are consistent with the results of Refs. [16-18]. The solid lines are the predictions of the full 



calculation including the anomalous terms. It is evident that the anomalous terms give non- 
trivial effects to the cross sections especially for the 7r + J/if> reaction. For comparison, we 
also calculate the normal processes considered in previous studies but with the additional 
diagrams with anomalous vertices. The results are given by the dot-dashed lines in Fig. 
|j. The difference between the dashed lines and dot-dashed lines are due to the anomalous 
processes, i.e., for n + J /if) — > D* + D,D + D* in the upper panel and for p + J/if> — > 
D + D,D* + D* in the lower panel. At low energy, i.e., \/s < 4 GeV which corresponds 
to E kin = s/s — Mj/^ — M w < 0.8 GeV, where Mj/^ and M w are the J/if> and pion masses 
respectively, the anomalous interactions reduce the total absorption cross section for n + J/if) 
by about 50 %. Our numerical results show that the n + J /if) absorption cross section is 
less than 10 mb in this energy region. However, the contribution of the anomalous terms 
to the p + J /if) cross sections is not so important except at the high energy region, which is 
thought to have many corrections from higher meson resonances exchanges. 

In Figs. [5] and [| we show the contributions from each diagram for each channel. For 
7r + J /if) — > D + D and ir + J /if) — > D* + D*, it can be seen that the contact diagrams 
(lc) and (3e) are important. However, considering the uncertainty related to the 4-point 
vertices, the results for this process are more qualitative. The more important and less 
uncertain process is the it + J /if) — > D* + D. The interesting observation here is that the role 
of the anomalous term [(2b) of Fig. |l| is quite non-trivial at low energies and even dominates 
at higher energies. The strong interference with the normal terms leads to the substantial 
difference in the cross section as shown in the upper panel of Fig. |j. But in contrast, as can 
be seen from Fig. || the role of the anomalous terms in the p + J /if) absorption is not so 
crucial except p+ J /if) — > D* + D,D + D* which is allowed only by including the anomalous 
vertices. Its effect becomes important at higher energies and is small at ^/s < 4 GeV. 
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C. Cross sections with form factors 



As pointed out by Ref. [|14||, the assumption of t-channel exchange of a heavy meson 



is hard to justify without including form factors. This is so because the range of heavy 
meson exchange is much smaller than the physical sizes of the initial hadrons. There are 
some approaches to include form factors in the J/ if> absorption processes and different form 
factors and cutoff parameters are employed |T7| , [T8|j36|| , which, of course, cannot be justified 
a priori. Microscopic models such as quark potential models may give us a guide for the 



form factors. Here, we employ the simple form of the form factors employed in Ref. |T7 
which are 



A 2 „ A 2 A 2 



A 2 + r 2 * A 2 + f 2 A 2 + f 



(32) 



where r 2 = (p x - p 3 ) 2 or (p 2 - p 3 ) 2 and f = [(pi - p 3 ) 2 + (p 2 - p 3 ) 2 }/2. F 3 is the form 
factor for the 3-point vertices and F4 for the 4-point vertices. The cutoff parameters should 
be chosen experimentally and may take different values for different vertices. However, 
because of the paucity of experimental information, we use the same value for all cutoff 
parameters as in the literature and investigate the dependence of the cross sections on the 
cutoff parameters. 

In Fig. we show a similar figure as in Fig. f| but with the form factors of Eq. fl32f ) with 
A = 2 GeV (left panel) and 1 GeV (right panel). It can be seen that the role of the form 
factors are more important at higher energy region and suppresses the overall cross sections. 
But the tendency shown in Fig. [| is still valid, i.e., the anomalous terms reduce the cross 
sections for 7r + J/if> at low energies through interference with the normal terms. Finally, 
Fig. [7] then shows that the cross sections for 7r + J /if; absorption process is about 2 ~ 6 mb 
depending on the cutoff parameter at \/s < 4 GeV. The p + J /if) cross section is very small 
at y/s < 4 GeV and has maximum value of 3 ~ 9 mb for intermediate energies up to 5 GeV. 



IV. SUMMARY AND DISCUSSIONS 

Knowing the magnitude and energy dependence of the J /if) absorption cross sections 
by hadrons is crucial to estimate the effect of J /if) suppression in RHIC due to hadronic 
processes. We have re-examined the J/ if) absorption by 7r and p mesons in a meson exchange 
model. The absorption processes are dominated by the D and the D* exchanges. By 
including anomalous parity interactions we found that additional channels and absorption 
mechanisms are allowed. Their contribution was found to be quite non-trivial especially 
for 7r + J /if) cross section. The calculated cross section for tt + J /if) at low energies was 
found to be only about one half of the previous calculations fll6Hl8f. With conventional, 
but arbitrarily chosen, form factors and cutoff parameters, the 7r + J/if> cross section was 
estimated to be 2 ~ 6 mb for \/s < 4 GeV. In contrast, the effect of anomalous terms was 
found to be weak in p + J/if> cross section at low energies. 

In addition to the uncertainties related to the determination of the coupling constants 
and form factors, our model calculation may be improved further. One such improvement 
would be the inclusion of the axial- vector Z) 1 (2420) meson exchanges, as mentioned in Ref. 
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T6| . Also at low energy, final state interactions are expected to improve the calculation. It 



would be also interesting to investigate other possible absorption channels as in Ref . [18| . For 
example, one may consider OZI-evading vertices such as the J/ip-p-ir coupling, which can 
allow new absorption channels with much smaller thresholds than the processes considered 
in this work, although such processes are expected to be suppressed. 
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APPENDIX: 



In this Appendix we briefly explain the way to obtain the effective Lagrangian (|l|) and @. 
As in the literature JT5|-|TS|] we start with the SU(4) symmetric chiral Yang-Mills Lagrangian 



27] 



£ X ym = ^ Tr {D,UD^) + 3- Tr [M (U + - 2)] 

o o 

- I Tr (F^F L ^ + F R F R ^) + m\ {A$A L * + A R A R ») , (Al) 

where U = exp(i2(p/ F n ) with the pion decay constant F w (« 132 MeV) and (p = TT a X a /\^2. 
The covariant derivative is defined as 

D,U = d,U - igAft + igUA% (A2) 

with left- and right-handed spin-1 fields A^ ,R , which are related to the vector and axial- vector 
mesons (= V£\ a />/$) and A, (= A a ^X a /^/2) by 

Al, = \iY, + A,) , A R , = I (V, - A,) . (A3) 
The field strength tensors are 

F L ,f = d,A L v ' R - d v A L * - ig [A^ R , A L ^ R ]. (A4) 



In obtaining the Lagrangian (|A1| ) we do not take any non-minimal terms into account for 
simplicity. 

It is straightforward to obtain the interaction Lagrangian for (f> and V, which leads to 



Ll p = -| Tr {V, (0^0 - ^#)} + | Tr (d,V u - d v V,)[V^ V] 

- ^ Tr [V„ 0] [V», <j>] + £ Tr [V„ V u }\ (A5) 
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where we have ignored the A-(j) mixing effects as they are of order of 1/M V . The above 
Lagrangian was obtained by keeping the axial- vector field. As in Ref. |[28| , one may gauge- 
out the axial- vector field. It is found that including proper non-minimal terms gives the same 
form for the effective Lagrangian as Eq. (|A5|) with different coupling constant |18], i.e., by 



substituting 2g for the coupling constant of Ref. |18| one can get the effective Lagrangian in 
the form of Eq. (|A5|) . 

The anomalous parity terms come from nine gauged Wess-Zumino terms, which read 



A, 



5Ci Tr [A L L 3 + A R R 3 ] - 5C Tr [(dA L A L + A L dA L )L + (dA R A R + A R dA R )R] 
+ 5CTr [dA L dUA R U* - dA R dU ] A L U] + 5CTr [A R U^A L UR 2 - A L UA R U ] L 2 ] 

+ T Tr[(A L Lr-(A R Rr] +5 CiTr[AlL + m 

+ 5Ci Tr [(dA R A R + A R dA R ) U^A L U - (dA L A L + A L dA L ) UA R U ] ] 
+ hCi Tr [A L UA R U ] A L L + A R U ] A L UA R R] 



+ 5CTr 



A%U*A L U - A\UA R U + -{UA R U^A L f 



(A6) 



where C = — iN c / 2A0n 2 with the number of color N c and L = dUU\ R = U^dU. We have 
used 1-form notation in writing the above Lagrangian. 

For simplicity we again ignore the A-<p mixing and use the Bardeen subtracted form for 
the anomalous terms. This gives 



C 



VP 
an 



' /X Tr(dVdV<f>)-^Tr {V{d<pf} 



ig 3 N c 



+ 



32n 2 F n 



Tr (y : 



+ 



Q n 2 F 3 

ig 3 N c 

32tt 2 F 7T 



Tr (VdVV(p) . 



(A7) 



Writing the pseudoscalar and vector fields explicitly and assuming pure cc for the J /if) 
wave function, we can obtain the effective Lagrangian (|J) and (0) as well as the SU(4) 
symmetry relations used in Sec. III. Note that by eliminating the axial-vector fields, the 
coupling constant of Vcfxfxp interaction becomes |28| 



My 



9v<t>4> 



Fl 



Assuming that the proper symmetry breaking terms make My 
F w — > F D in the above relation, we get g^ DD7V w 18.0 GeV- 3 , 
listed in Table |. 



(A8) 

-> M D *, g v<M> -> g D , D ^ and 
which is close to the value 
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TABLES 



coupling constant method value 

g D . Dn D* -> Dtt and QCD sum rule 8.84 || 

9ipDD 1 DD coupling and VDM 7.71 g| 

g DDp jDD coupling and VDM 2.52 |f 

g^ D , D , -fD*D* coupling and VDM 7.71 ||Jl 

g D . D . p -fD*D* coupling and VDM 2.52 |Jjl 

9ipD*Dw SU(4) symmetry relation 33.92 |TJ] 

9ipDDp SU(4) symmetry relation 38.86 [17] 

SU(4) symmetry relation 19.43 [17] 



9d*D*tt 


heavy quark spin symmetry 


9.08 GeV" 1 




D* — ► D7, VDM, and quark model 


8.64 GeV" 1 


9d*d p 


D* — * Dj, VDM, and quark model 


2.82 GeV" 1 


9^DDtt 


SU(4) symmetry relation 


16.00 GeV~ 3 


9ipD*D*n 


SU(4) symmetry relation 


38.19 GeV- 1 




SU(4) symmetry relation 


38.19 GeV- 1 


9ipD*Dp 


SU(4) symmetry relation 


21.77 GeV" 1 


h J ij ) D*Dp 


SU(4) symmetry relation 


21.77 GeV" 1 



TABLE I. Determination of the coupling constants. 
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FIGURES 




(la) (lb) 



D D* D d" 



Sly " Jlf 
(2a) (2b) 



D D D D D D D 





(3a) 



(3b) 



(3c) (3d) (3e) 

FIG. 1. Feynman diagrams for J/ip absorption processes by pion: (1) J/ip + ir — > D + D, (2) 
J/i/> + tt -> D* + D, and (3) J/i/> + vr -> D* + D*. The process J/i/> + vr -► L> + 5* has the same 
cross section as (2). 




J/y p J/y p J/y p JA|r p J/y 

(4a) (4b) (4c) (4d) (4e) 



D D* D D* D D* D D* 




p J/y p J/y p J/y p J/y P J/V 

(5a) (5b) (5c) (5d) (5e) 



D D D D 




(6a) 



J/y p p 
(6b) 



JAP p 



(6c) 




FIG. 2. Feynman diagrams for J/ip absorption processes by p: (4) J/tp + p — > D + L>, (5) 
J/V> + p -»■ -D* + A and (6) J ftp + p ^ D* + D*. The process J/^ + p -> Z) + £>* has the same 
cross section as (5). 
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io 3 i | 




(GeV) 

FIG. 3. Cross sections for ir + J/ip and p + J/tp. In the upper panel, the dotted, dashed, and 
dot-dashed lines correspond to Diagrams (1), (2), and (3), respectively. In the lower panel, the 
dotted, dashed, and dot-dashed lines correspond to Diagrams (4), (5), and (6), respectively. The 
solid lines are the sum of all processes. 




Vs (GeV) 

FIG. 4. Cross sections for ir + J/ip and p + J/tp- The solid lines are our model predictions and 
the dashed lines are those obtained without the anomalous terms. The dot-dashed lines correspond 
to 7r + J /if) — ► D* + D, D + D* in the upper panel and p + J/tp -> D + D, D* + D* in the lower 
panel in full calculation. 
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Jt+JA|/-»D+D JT+JA|f->D*+D,D+D* Jl+J/\|f->D*+D* 



io J 



10" 



10 1 

t 10° 
D 



10 



10" 



10 



-3 




(la,lb) 

(ic) m 



sum 

i 1 1 1 i I i i i i I i i i i 



I I I I | I I I I | I I I I: 




(2a) 

---- (2b) 
--- (2c) 
■ — (2d) ^ 

sum 

I i i i i I i i i i 



I I I I | I I I I | I I I I 




(3a,3b) ; 

(3c,3d) _ 

--■ (3e) 

sum 

i i i i I i i i i 



3.5 4 4.5 3.5 4 4.5 3.5 4 4.5 5 



Vs (GeV) Vs (GeV) Vs (GeV) 

FIG. 5. Contributions from each diagram to the cross sections for it + J/ip. 



p+JA|/->D+D p+JA|/->D*+D,D+D* p+JA|/-»D*+D* 



10 M i | i i i i | 



io 1 



t 10° 



10" 



10 



-2 



m — i i i i i — i i =|= i i i i i — r-r-| — i i i i =|= i i | i i i i | — n - ! — r 



(4a,4b) 
-- (4c,4d) =E 

■- (4e) X 
— sum 
i i I i i i i I i i i i 




i i i i l i i i i 



I 1 1 1 1 I 



(6a,6b) 
(6c,6d) E 
(6e) 
sum 



i i i i I i i i i 



4 4.5 5 4 4.5 5 4 4.5 5 
Vs (GeV) Vs (GeV) Vs (GeV) 



FIG. 6. Contributions from each diagram to the cross sections for p + J/tp. 
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A = 2.0GeV A=1.0GeV 




Vs (GeV) Vs (GeV) 

FIG. 7. Cross sections for ir + J/ip and p + J/tp with A = 2 GeV (left panel) and 1 GeV (right 
panel). The notations are the same as in Fig. ||. 
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